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ABSTRACT 


An  equilibrium  slab  model  of  a  barium  cloud  is 
derived.  This  model  has  finite  length  in  the  mag¬ 
netic  field  direction  and  allows  one  to  study  the 
stability  of  configurations  with  various  degrees 
of  ionospheric  short-circuiting.  It  reduces  to 
previous  models  for  the  case  of  no  short-circuiting 
but  is  quite  different  in  the  more  general  case. 

The  combined  influence  of  finite  Larmor  radius 
corrections  and  inertia  on  the  stability  of  the  in¬ 
finite  slab  model  is  studied.  A  simple  dispersion 


relation  is  obtained. 
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INTRODUCTION 


In  previous  semi-annual  technical  reports,  we  have 
developed  a  three-dimensional  model  of  the  equilibrium 
(i.e.  pre-striation)  growth  of  a  barium  cloud  in  the  pre¬ 
sence  of  a  conducting  background  ionosphere.  Early  stu¬ 
dies  considered  the  case  of  a  highly  conducting  back¬ 

ground  with  complete  short-circuiting.  In  this  case  there 
is  no  polarization  of  the  cloud  and  a  diffusing  elliptical 
gaussian  density  results.  Later  work(3»4)  extended  these 
results  to  include  finite  background  conductivity  with 
polarization  effects  and  steepening  of  the  backside. 

Stability  calculations  (5)  have  been  carried  out  (in  the 
case  of  the  highly  conducting  background)  by  studying  the 
growth  of  small  perturbations  about  the  equilibrium  de¬ 
veloped  above.  These  numerical  calculations (^)  have  shown 
that  the  fastest  growing  eigenraode  is  concentrated  on  the 
rear  side  of  the  cloud.  Calculational  difficulties  have 
not  yet  enabled  one  to  establish  the  stability  threshold 
(or  whether  one  actually  exists)  nor  has  it  yet  yielded  a 
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this  point  is  in  progress. 
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Meanwhile,  some  information  can  be  gleamed  from  slab 
models,  although  there  are  troublesome  questions  concerning 
the  validity  of  application  to  real  three-dimensional  clouds. 
The  E  X  B  instability  was  first  derived  in  1963  using  such 
a  model  by  Simon(^)  and  by  Hoh(^)  with  parallel  to  the 
density  gradient.  This  analysis  yielded  a  threshold  value 
of  £?  and  a  critical  wavelength  which  was  the  first  mode  to 
go  unstable  (and  the  fastest  growing  at  any  later  "linear" 
stage).  Linson  and  Workman  extended  this  analysis  to  in¬ 
clude  a  component  of  the  field  perpendicular  to  the  density 
gradient.  In  their  treatment,  they  take  the  limit  of  col¬ 
lisionless  electrons  £b  initio  with  the  result  that  there 
is  no  threshold  for  instability  onset  and  all  wavelengths 
grow  at  about  the  same  rate.  This  result  is  disturbing 
because  it  seems  at  variance  with  observed  time-delay  for 
striation  onset  and  because  the  analysis  gives  no  guide 
for  predicting  the  size  of  observed  striations. 

There  are  a  number  of  effects  which  could  restore  a 
threshold  to  the  slab  models  above.  These  include  effects 
of  ion  inertia  and  finite  larmor  radius  which  are  discussed 
in  Sect.  Ill  below.  Another  effect  is  the  internal  equi¬ 
librium  electric  field  that  developes  in  the  plasma.  In 
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the  infinite  slab  model,  such  as  that  of  Linson  and  Workman, 
this  necessarily  points  towards  the  center  of  the  slab  in 
order  to  keep  ions  from  diffusing  out  faster  than  electrons. 
This  field  direction  is  destabilizing  on  the  basis  of  the 
Simon-Hoh  model.  In  Sect.  II,  we  develop  a  model  of  a  slab 
with  finite  extent  along  the  magnetic  field  and  allow  for 
a  background  of  arbitrary  conductivity.  The  resulting  field 
can  either  point  in  or  point  out  depending  on  the  degree  of 
short-circuiting  by  the  background.  A  stability  analysis  of 
this  more  general  slab  model  is  now  in  progress. 
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II.  SLAB  MODEL  IN  A  FINITE  CONDUCTING  BACKGROUND 


The  model  we  consider  is  sketcltf.c  ia  Figure  1.  The 
cloud  is  represented  by  a  slab  of  finite  extent  along  the 
magnetic  field.  We  shall  assume  no  variation  of  the  equi¬ 
librium  quantities  in  the  ^  direction,  \J  X  &  z  0  then 

implies  const,  *  Eof 

As  shown  in  a  previous  report (3),  if  we  expand  the 
current  divergence  equation 


V  ■  7 


(1) 


in  the  small  parameter  Dtl  / T)„  we  obtain  in  the  zeroth 

order  approximation,  the  electric  field 


-A^  -  vfi.ro 


(2) 


where  M  is  the  equilibrium  density. 

To  determine  fyOj.)  we  shell  integrate  eqn.  (1)  along 
a  magnetic  field  line  from  one  boundary  of  the  ionosphere 
to  the  other  boundary.  Using  the  condition  that  the  current 
vanishes  at  the  boundaries,  we  have 


— .  0 


(3) 


In  the  large  (j. cT  )_  limit 


i  = 


DIvJn  ^  (tslBj.) -/jUOQ  6.  v-V.e 
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where  the  symbols  all  have  thr  some  meaning  ss  in  previous 
reports. 

If  we  substitute  eqn.  (2)  into  eqn.  (4)  and  integrate 
the  equation  along  a  field  line,  we  find 


Zffc 

.?x 


_l<p£Z 

e  ^ a 


o 


(5) 


In,  arriving  at  eqn.  (5)  we  set  k!  ■=.  fj, -(■  where  NA  is  the 


A 


ambient  density,  and  have  defined 

/  Jf/i 
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Nc  d } 
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S.  =  e  )  A <1  N 'A  d j 

-L 

We  also  assume  that  (  ClT is  independent  of  J.  within  the 
length  of  the  cloud. 

Integrating  eqn.  (5),  we  find 

-  (  t>Uaj,C  )  &  -  (pi  %  +  )  41 


—yih  (dlT),  ?dc  =•  tc-M*. 

To  evaluate  the  constant*  we  consider  a  field  line  far  away 


from  the  cloud,  Where  hj '  z  0  ,  and 


Jit 

}  x 


f.c X 


we  find 


/ 


2  it 


/!  ‘-CX 


const. 


We  have  then 


n  ^ 

?x 


o(  —  (.1  - ct ) 

^  )jj.  /  7A/6ax  **  (nz^ 


where 


c*w  =  Zfi/ie/Ait 3C&)+  tfA  1 


(6) 


(7) 


is  the  ratio  of  Pederson  conductivity  integrated  along  a 
magnetic  field  line  not  passing  through  the  cloud  to  that 
integrated  along  a  field  line  through  the  cloud. 

It  follows  from  eqns.  (2)  and  (6)  that 


+  (Xtox  -r  ( / '  & )  .OtL 

(.1X1 


(8) 


If  the  background  ionosphere  is  assumed  to  be  infinitely 
conducting,  then  £X-?  j  and  we  have 


r*  (0  *  K  I  <)  /y  I  {Z 

£  *  =  -  4  ;<j  7X  L^' 


in  agreement  with  the  electric  field  used  in  earlier  reports 
(cf.  eqn.  (9)  of  Reft  3). 

If  the  cloud’s  contribution  dominates  in  the  expression 
for  the  integrated  conductivity,  then  (X ->  C .  As  the  back¬ 


ground  density  is  negligible  we  have  ^ 


i  >  kJ< 


and  then 
£ 


a  > 


jzL  i-  id  e 

jui  *  (a-o+ 
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This  is  the  electric  field  used  in  the  slab  model  of  Linson 
and  Workman (8). 


The  expression  for  the  electric  field  to  the  next  order 
in  the  expansion  in  the  small  parameter  DJj / D^\  is  ob¬ 
tained  by  iteration  from  eqns.  (1)  and  (2).  We  find 

-  C kU^)'L£ 7  #  (9) 

where  £  -  (vj - Vj  *  MjuUjJI  ?l/(Jbli+)A7j) 

K  -=■  ( M  j  ) ;/(.  JUti  -V-yU  U  ) 

L  =  f tav*/Aj  -  J / ; 


To  find  the  equilibrium  density,  we  substitute  eqn.  (2) 
into  the  ion  continuity  eqn.  and  obtain 

-  (Vl-tA  /-U  )  -(Tj^ALlIP 

Making  use  of  eqn.  (6)  ws  have 

#  -  ( 0// ♦  v* ^  ^  - 1 « 


(10) 


'5  '  ** 


(11) 


If  the  ambient  density  is  negligible  inside  the  cloud,  then 

/Vii  A/t .  We  obtain  the  solution  for  N/c 
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(12) 


S L 


H:  rrrzTTf  ^  $  /^n  ^  ) 


<£)//  1=5  s  A /J  ii 

and  hjc  is  a  solution  of 

4>c  "  (  ^1-^A/Jl  ) 


where 


(13) 


In  the  case  of  an  infinitely  conducting  background, 
CX^  /  ,  we  have 


^  -  lvPajuI)  %  ^yi  2$  f7 


^K'c 

~Tx 


0 


The  solution  for  an  instantaneous  line  source  is 

S  O 


^  Tnzkrl  ex/>Hx-*0%^0 


M  411  by i 


(14) 


where 


£)x~viJ  Ay i 

vt  =  J 

We  see  that  the  strong  short-circuiting  will  allow  the  cloud 
to  diffuse  at  the  ion  rate  and  the  cloud  as  a  whole  moves 
with  the  ion  drift  velocity  in  the  ambient  electric  field. 

In  the  case  of  dominant  cloud  conductivitv 


Eqn.  (13)  becomes 


7>Xlc  ,  Cf  3/vt  _ 

7-1  -7x  - 


0 
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and  then 


(15) 


a£  =.  a/c  ( 't'U't) 

The  cloud  as  a  whole  moves  wi  th  the  velocity  U0 

Had  we  included  electron  collisions,  the  density  profile 

would  also  diffuse  at  the  electron  rate. 


III.  ION  INERTIA  AND  FINITE  LARMOR  RADIUS  EFFECT 

The  results  in  both  references  1  and  8  show  that  there 
is  no  threshold  for  instability  in  the  limit  of  collisionles 
electrons  and  with  no  variation  of  the  perturbation  in  the 
magnetic  field  direction.  This  cannot  be  generally  correct 
in  the  limit  of  collisionless  ions  as  well,  since  in  that 
case  one  can  remove  the  constant  external  electric  field 
by  a  uniform  frame  change  with  no  alteration  of  the  plasma 
dynamics.  It  is  well  known  that  the  inclusion  of  previously 
omitted  inertial  terms  in  the  equations  of  motion  is  one 
resolution  of  this  paradox.  However,  the  quantitative 
threshold  for  instability  onset  is  too  low  to  account  for 
cloud  observations. 

Another  mechanism  with  a  similar  effect  is  inclusion 
of  finite  Larmor  radius  correction  in  the  fluid  equations 
of  motion.  For  completeness,  and  because  some  of  the  tech¬ 
niques  are  unfamiliar,  we  have  calculated  the  linear  dis¬ 
persion  relation  with  both  included. 

Roberts  and  Taylor (9)  have  shown  that  one  can  include 

the  finite  Larmor  radius  effect  by  proper  modification  of 

the  ion  pressure  tensor  in  the  magnetohydrodynamic  equation. 

The  ion  momentum  equation  including  inertia  and  finite  Lar- 
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mor  radius  takes  the  form 


</t/* 

where 


-r  f'  \/  /  &  Y)  A'l  v  ^  I  \ 

=:  -KTvn  -#•  +necL-^-~ 

~  c  ~c 


(16) 


v=^a2a 


3V^ 


^  J 


where  A  is  the  ion  Larmor  radius. 

The  ratio  of  the  finite  Larmor  radius  term  to  the 
collisions!  term  is  Oxlcd/ji  )  where  J  is  some  characteristic 
length.  We  shall  limit  our  discussion  to  a  slab  model  with 
no  short-circuiting  (#? (?)  and  uniform  in  the  j  direction. 

The  equilibrium  ion  velocities  are 

v/=  -  vlfr# 


C 

0 


L\ 


V*  o 

where  we  have  inserted  the  equilibrium  electric  field  ob¬ 
tained  in  the  previous  section.  Note  that  the  equilibrium 
velocities  are  constant  in  time  and  uniform  in  space.  It 
is,  therefore,  consistent  to  omit  the  inertial  and  the  finite 


Larmor  radius  terms  in  the  equilibrium  momentum  equation. 
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The  linearized  momentum  equation,  in  the  frame  moving 
with  velocity  UD  ,  and  with  1 s  ' 


-  i w  lk i  -  ik  p,  <t’+£c £*&- a?) 

Here  4 

^ f fox ) V/ x  —  L^$'d 

In  arriving  at  eqn.  (17),  we  assumed  that  *  '7 v 

is  constant  and  Fourier  analyzed  the  perturbed  quantities. 

-y* 

It  is  not  difficult  to  solve  for  V/x-  and  V,y  in  terms 
of  0  and  j  and  substitute  the  expressions  into  the 
linearized  ion  continuity  equation  in  the  primed  coordinates. 
We  find 

{-  Uu-r  L  iV  Jjj  t  i  ' 

+Mf^$*##r***«|*#* . .  „„ 


here 


p=  i ^ ) ] 


The  inertia  term  and  finite  Lormor  radius  effect  is  negligible 
in  the  electron  equation  of  motion.  The  linearized  electron 


equation  is  easily  seen  to  be 

!7i?  I?  ■+  1 V*"  t  ? 


l-i*  '  ikM 44  }-/ 


-f  J-  £ 

t 


0  ^  C 


(19) 
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From  eqns.  (18)  and  (19)  we  obtain  the  dispersion  relation 


In  the  limit  of  vanishing  Larmor  radius  ft-o  (hence  l>zo) 
and  with  the  assumption  tojlX^I*  eqn.  (20)  becomes: 


Therefore 

«  n-r  A,  -+  AiriAi/— 

0.!^  <y.  — -  -fi;ax -n^rnt  (21) 

In  this  limit  we  have  the  result  of  Linson  and  Workman^) 
when  proper  Doppler  shift  in  a)  is  allowed  for  between  the 
moving  frame  and  the  neutral  rest  frame. 

If  we  take  the  collisionless  limit  of  eqn,  (20)  we  find 

—  IujL  '-h1"  ^r)\  =  o 


Note  that  the  dependence  on  U0 
small  terras  of  order  (a/^)  and 
equation  divided  by  reduces 


vanishes. 


✓  /->  / 1\ 
\wa  / 


to 


If  we  neglect 

1^  4  A  i«  l*  A  ^ 

ui^uci.  ,  uic  cauwc 
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0  (22) 

There  are  two  low  frequency  inodes  (co^fl)»  One  has  ojz  0  and 
the  other  is  a  purely  oscillating  mode  with  U)  &  2^/d 
The  third  root  has  a  high  frequency  of  order  Ji  J Cl  and  is  not 
relevant  to  the  present  discussion. 

If  we  define  the  equivalent  gravitational  force  -t Uo/'T 
and  then  take  the  collisionless  limit  of  eqn.  (20),  neglecting 
small  terms  of  order  ( c\-A)1 ,  ( Q./d f  and  (■ o/Ci  we  have  the  dis¬ 
persion  relation  for  waves  propagation  perpendicular  to  the 
density  gradient  &) 

U)1-  UPAyrf-  =  0  (23) 

This  also  agrees  with  that  obtained  by  Roberts  and  Taylor^) 
when  the  difference  in  the  sign  of  u)  is  allowed  for.  Thres¬ 
hold  estimates  are  now  being  investigated. 
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FIG.  I  FINITE  SLAB  IN  A  MAGNETIC  FIELD 


